global QV qvout temp;

qvx = [
0
15.9
31.7
47.6
63.5
79.3
95.2
111
127
143
159
175
190
206
222
238
254
270
286
302
317
333
349
365
];

qvy = [
0
0.032
0.03
0.035
0.03
0.028
0.008
0.015
0.015
0.008
0.008
0
0.01
0.01
0.01
0.01
0.01
0.015
0.02
0.025
0.025
0.02
0.025
0.03
];

qvt = 0:1:365;
qvt = qvt';
QV = interp1(qvx,qvy,qvt,'linear');


qvoutx=[
0
36.5
73
110
146
182
219
256
292
328
365    
];

qvouty=[
0
0.032
0.03
0.032
0.02
0.015
0
0
0.015
0.025
0.025    
];

qvoutt = 0:1:365;
qvoutt = qvoutt';
qvout = interp1(qvoutx,qvouty,qvoutt,'linear');




tempx=[
0
30.4
60.8
91.2
122
152
182
213
243
274
304
335
365
];

tempy=[
0
0
2.2
5.8
12.7
16.5
19.8
22.4
18.7
15.2
9.8
5.4
1
];

tempt = 0:1:365;
tempt = tempt';
temp = interp1(tempx,tempy,tempt,'linear');


clear qvx qvy qvt qvoutx qvouty qvoutt tempx tempy tempt











